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INTRODUCTORY  REMARKS 


The  material  given  in  Technical  Report  No.  15(6l)  represented  a  "Long 
Abstract"  as  required  by  ACS  regulation  of  papers  to  be  presented  at  a 
national  meeting.  Jbllovlng  suggestions  made  after  the  meeting  at  Miami 
Beach  where  Dr.  Yajnik  presented  the  material,  it  was  decided  to  make  the 
material  available  to  the  general  public  by  publication.  This  decision, 
on  the  other  hand,  entailed  a  rewriting  of  the  paper,  an  appreciable  increase 
in  documentary  material  beyond  that  originally  available  and  additional  j 
checks  of  the  computations  made ,  The  resulting  considerably  revised  and 
enlarged  edition  which  is  now  being  submitted  for  publication  in  the  J. 
Polymer  Sci., differed  sufficiently  from  Technical  Report  No.  15(6l)  that 
its  re-issuance  as  Technical  Report  No.  15A(6lA)  scored  to  be  warranted. 


ELECTRONIC  CO-IPUIATIONS  OP  LIGHT  SCATTERING  JUNCTIONS 
FOR  HETJRODISPERSE  SYOTJMS  OF  ISOTROPIC  SPHERES 


M.  Yajnik,  J.  Witeczek  and  W.  He  Her 
Chemistry  Department 
Wayne  State  University 
Detroit,  Mich.  48202 


I.  INTRODUCTION 

An  important  problem  in  the  characterization  of  disperse  systems,  such 
as  polymer  latic*****,  emulsions,  etc.,  ia  the  determination  of  their  particle 
size  distribution.  Various  experimental  techniques  have  been  used  for  this 
purpose.  Light  scattering  methods  are  among  the  most  promising  ones.  They 
have  the  advantage  over  other  methods  (e.g.  centrifugation,  electron  microscopy) 
in  that  they  are  relatively  simple  experimentally,  do  not  JTect  the  state  of  the 
systems  and  are  very  sensitive  to  rapid  changes  in  particle  size.  It  is  therefore 
possible,  in  principle,  to  obtain  "instantaneous"  size  distribution;  and  to 
study  their  variation  with  time.  The  advent  of  computers  has  greatly  enhanced 
the  attractiveness  of  lighr.  scattering  methods  for  the  determination  of  size 
distribution  curves  since  thereby  the  formerly  very  appreciable  amount  of  time 
required  for  calculations  has  been  reduced  drastically. 

The  purpose  of  this  paper  is  to  show  how  the  use  of  a  computer  can  alloy 
one  to  obtain  the  size  distribution  of  relatively  large  spherical  particles 
from  light  scattering  data.  We  will  be  concerned  specifically  with  the 
application  of  the  Mia  theory  because  the  systems  investigated  contained  isotropic 
spherical  particles  whose  refractive  index  differs  appreciably  from  that  of 
the  medium.  The  computer  was  used  for  the  solution  of  two  problems;  firstly, 
e. .-pensive  tables  of  light  scattering  functions  were  prepared  for  a  given  type  of 
distribution  to  he  discussed  below.  Secondly,  after  the  experimental  light 
scattering  data  had  been  obtained,  the  computer  was  used  in  order  to  find  the 
best  fit  between  the  experimental  data  and  the  tabulated  the  or  <tical  values . 

Only  the  first  part  of  the  computer  application  will  be  considered  here. 

II.  DEFINITIONS  AND  HISTORICAL  REMARKS 

While  it  is,  in  principle,  possible  to  compute  particle  size  distributions 


without  Introducing  any  assumption  as  to  the  type  of  distribution  function, 
this  approach  Is  at  present  not  practical.  Therefore,  it  has  been  assumed1 
that  the  distribution  Is,  In  agreement  with  most  distributions  in  latices 
and  emulsions,  unimodal  and  positively  shewed.  The  function  Is 


f(r)  =  r  >r. 


=  0 


T*$r, 


Here  cf(r)dr  is  the  number  of  particles  per  unit  volume  with  radii  between  r 

and  r  +  dr,  c  is  a  normalisation  factor;  and  s  and  rQ  are  the  distribution 

parameters:  s  is  proportional  to  the  distribution  width  and  r  is  the  radius 

o 

of  the  smallest  particles  n resent  in  statistically  significant  numbers.  It  is 
advantageous  to  introduce  .  dimensionless  parameters 
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where  X  iB  the  wavelength  of  the  incident  light  in  the  medium. 

There  are  several  light  scattering  functions  which  can  be  used  in  the 
evaluation  of  particle  size  distributions. 

The  two  methods  developed  previously  in  this  laboratory  involved  the 
determination  of  the  spectra  of  turbidity  and  of  the  spectra  of  the 
scattering  ratio  at  an  angle  of  observation  of  90° •  The  scattering  ratio 
«•-  V1^  ,  vhocc  lf|  and  1^  are  the  total  intensities  of  light  scattered  from 
an  incident  linearly  polarized  beam  of  unit  intensity  with  the  electric 
vector  vibrating  in  the  plane  of  observation  and  perpendicularly  to  the 

plane  of  observation  respectively.  We  shall  consider  here  the  scattering 

c 

ratio  for  a  series  of  angles  of  observatio  in:luc,J:ix;  .O  .  Cotuequently, 

the  variation  of  5*  with  the  angle  of  observation  can  be  used  as  the  third 

* 

criterion  in  the  determination  of  particle  size  distributions  . 

Finally,  as  a  fourth  potentially  useful  argument,  the  angular  variation 
of  In  and  I_k  ,  themselves,  shall  be  considered. 


•Closely  related  to  this  is  the  use  of  the  angular  variation  of  the 
polarization  ratio  evaluated  by  Ker&er  et  al*  for  certain  relatively 
large  refractive  index  ratios  applicable  to  aerosols. 
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Considering  a  system  containing  optically  homogeneous,  non-absorbing 

spheres  with  a  distribution  in  size  f(r),  the  quantities  enumerated  are, 

3 

on  using  the  Mie -equations  , defined  as  follows: 

I,  -  <30 

lxm  ^  f  L(*.!f>»-)'f(r)dr  (3b) 

and 

5-  =  I  (30 

Jr  (“.  Ji «-)  •  f(r)  ctr 

Making  U3e  of  equations  (2)  we  get  the  following  expression  for 
the  scattering  ratio: 


_  j  1  (4) 

Ihe  quantities  i„  (a  ,  7 ,  m)  and  ij_  (a  ,  7 ,  m)  in  equations  (3) 
and  (4)  are  the  basic  Mie  functions  for  angular  scattering.  They  depend 
on  a,  on  the  refractive  index  of  the  spheres  relative  to  that  of  the 
medium,  m,  and  on  the  angle  of  observation  7 ,  measured  with  respect  to 
the  reverse  direction  of  the  incident  beam.  They  are  defined  as  follows 


.0,1  =  |  J  [A„r,,(««|f)  +  (<*)[)]  |  <50 

^  =  j  ZI  [^K(u«^)+  B„Tj«s!j)]|  '50 
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where 


A„=H) 


A, 


2  n  +  l 
a  (w  +  I  ) 


£,fr)S.7g)-m,S»gJ|i) 


2n  +  I 

■'(”<■')  0;  («)S„  (ft)-  m  $4, (*)$.'&) 


ft  = 

S  are  the  Rlccati-Bessel  functions,  ©  are  complex  functions  related 
n  '  n 

to  the  Hankcl  functions.  'TT  and  T*  are  related  to  the  associated  Legeudre 

n  n 

polynomials.  Puller  description  of  the  quantities  which  enter  in  the  Mie 

4 

equations  has  been  given  elsewhere  . 

The  first  systematic  attempt  to  compute  solutions  of  the  Mie  equations 
(5a)  and  (5b)  for  specific  values  of  the  parameters  a,  y  and  m  and  for 

5 

aonodisperse  systems  by  Blumer,  in  1925  had  to  make  use  of  manual  computations. 

The  most  recent  results  obtained  by  manual  computations  dealing  r.,rain  with 

monodlsperse  systems,  were  those  by  Pangonis.^  They  were  calculated  with  a 

conventional  de6k  calculator.  Since  it  is  necessary  to  take  at  least  a  terns 

and,  preferably,  (1.43  a  +  5)  terms  into  the  summations  defined  in  equations 

(5)  the  practically  possible  scope  of  computations  was  severely  limited^ 

however,  unless  sl  uv  .•  >.i:.Lly  large  amount  of  tine  was  to  be  devoted  to  the 

task.  Thus,  the  computations  by  Pangonis  referred  to  above,  required  the 

better  part  of  one  year.  Only  since  the  advent  and  more  general  availability 

of  electronic  computers  has  it  become  possible  to  evaluate  the  Mie  scattering 

4  6  7 

functions  without  major  restrictions  in  parameter  values .  *  *  Tbe  need  for 

computations  by  an  electronic  computer  became  even  more  imperative  on  using 

the  Mie  theory  for  the  determination  of  distribution  functions.  Here  manual 

computations  are  virtually  impossible  because  of  the  prohibitive  time 

0 

required.  Stevenson  and  Heller  in  a  first  attempt  in  this  direction  have 
published  electronically  obtained  scattering  functions  for  heterodisperse 
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systems  for  S'- values  at  the  single  angle  7  = 

in.  THE  PROGRAMMING  TECHNIQUE  AND  PERFORMANCE  CHECKS 

The  program  as  developed  for  the  present  purpose  was  written  mostly 

in  Fortran  for  the  XEM  707^  computer  and  consists  of  several  parts.  They 

are  concerned  with  (l)  the  generation  of  Bessel  functions  and  Legendre 

polynomials  (2)  the  evaluation  of  Mie  functions  defined  by  equations  (5),  and 

(-3)  wit!  ‘Usir  integration  for  specific  p  and  q  values.  The  generation  of  the 

Riccuti-Beasel  functions  was  carried  out  by  downward  recursion  since  this 

o 

has  been  found  to  be  more  accurate  i'or  higher  orders.  In  thi3  connection  the 
method  based  on  the  expression  of  ratios  of  ife-.sel  functions  in  terms  of 
continued  fractions  has  been  used.10 

On  carrying  out  the  integration  for  a  given  pair  of  p  and  q  values,  it 
was  of  course,  necessary  to  cut  off  the  process  of  integration  at  that  a 
value  above  which  the  asymptotically  decreasing  contribution  of  higher  O- values 
to  the  total  scattering  is  negligible.  It  was  found  that  the  desired  accuracy 
of  the  scattering  foncttor  iigiuws)  was  not  affected  if  the  cut-off 

was  made  at  that  Ct-value  at  which  the  value  of  the  distribution  function  f  ,:*) 
has  decreased  to  less  than  10"'>. 

One  of  the  major  problems  to  be  resolved  was  the  choice  of  the  best 
method  of  integration.  For  this  purpose  extensive  experimentation  was 
undertaken  by  means  of  the  computer,  the  objective  being  to  obtain  a  sufficiently 
accurate  result  without  having  to  use  an  inordinate  amount  of  machine  time.  Two 
classes  of  methods  were  used:  (l)  The  Newton-Cotes  formulae,  and  (2)  The 
Gaussian  Quadrature.  Since  the  former  are  equal-interval  formulae  they  are 
simpler  to  use,  particularly  if  the  values  of  the  integrand  are  already 
available.  This  is  due  to  the  fact  that  one  set  of  data  can  be  used  for  al.1 
Integrations  (for  any  p  and  q  values).  For  preliminary  checks,  some  manual 
computations  were  performed  using  both  the  trapezoidal  and  Simpson's  rule. 

Both  performed  satisfactorily.  Simpson's  rule  was  considered  as  preferable 
however,  and  adopted  because  of  the  oscillatory  character  of  the  functions 
at  large  ex-values.11 

It  is  true  that,  in  principle,  the  Gaussian  Quadrature  should  give  more 
reliable  results  if  the  number  of  Intervals  between  the  integration  limits  of 


5 


the  independent  variable  is  made  sufficiently  large.  Therefore,  this 

alternate  method,  was  investigated  also  more  closely.  The  process  ox' 

Integration  was  allowed  to  continue  automatically  for  smaller  and  smaller 

rvi  ’1  Isions  until  the  difference  between  two  successive  results  was  less 

than  10-^.  The  general  procedure  is  schematically  shown  in  Figure  1.  which 

# 

was  prepared  with  a  program  using  the  autochart  programming  system. 
Unfortunately,  the  iterations  anr.  the  necessity  of  evaluating  the  Hie  functions 
in  each  instance  took  up  a  major  amount  of  computer  time.  Therefore,  this 
method  is  used  at  the  present  time  for  spot-checks  only. 

A  comparison  of  the  results  obtained  with  the  two  methods  of  integration 
is  given  in  Table  I.  where  the  bracketed  data  are  taken  from  Reference  8;  in 
view  of  the  fact  that  the  computing  techniques  improved  in  the  meantime,  the 
*C'  -i  •  h  'he  present  values  is  satisfactory.  It  is  apparent  that  a 
aecroico  in  the  jx-interval  fra.i  0.2  to  0.1  does  not  materially  change  the 
result  obtained  with  Simpson's  rule,  while  a  decrease  is  observed  in  the 
results  obtained  with  the  Gaussian  method  on  changing  from  5-node  to  9-node 
computations . 

A  comprehensive  survey  of  the  performance  of  the  Gaussian  Quadrature  is 
given  in  Figure  2  where  the  number  of  nodes  is  the  independent  variable  and 
the  number  of  panels  is  the  parameter. 


*  The  flow  chart  given  in  Figure  1  illustrates  in  principle  also  the  computing 
procedure  on  using  Simpson's  method  except  that  the  absence  of  nodes 
eliminates  the  respective  part  of  the  process.  In  addition,  changing  intervals 
between  successive  a  values  is  here  much  simpler. 

**  The  number  of  nodes  define,  the  fixed  number  of  intervals  into  which  each 
panel  is  subdivided.  A  panel,  on  the  other  hand,  defines  onr  of  the  N  parts 
Into  which  the  total  interval  of  integration  is  originally  divided.  After 
the  first  computations  the  number  of  panels  is  increased  automatically 
to  2N,  4N,  etc.,  the  number  of  nodes  being  in  each  instance  constant.  The 
machine  was  programmed  to  Increase  the  number  of  nodes  systematically  after 
each  such  operation. 
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It  is  apparent  that  it  is  necessary  at  a  given  p  and  q-value,  to  have  a 
product  of  at  least  60  for  the  number  of  panels  and  number  of  nodes  in  order 
to  obtain  a  result  which  does  not  change  further  on  increase  of  panels  or 
nodes.  It  is  new  Interesting  to  note  that  the  nearly  asymptotic  value 
obtained  on  using  a  sufficiently  large  number  for  the  product  of  nodes  and 
panels  is  close  to  the  result  obtained  on  using  Simpson’s  r jIc  provided  the 
number  of  intervals  in  the  latter  instance  is  close  to  100.  The  respective 
asymptotic  values  differ  by  somewhat  less  than  1.5$. 

The  satisfactory  convergence  of  the  results  obtained  by  the  two  methods 
(Simpson's  rule  and  Gaussian  Quadrature) on  meeting  the  proper  requirements 
with  regard  to  the  number  of  intervals,  zru.1  products  of  panels  and  nodes, 
respectively,  represented,  in  principle,  sufficient  proof  for  the  reliability 
of  the  results  obtained.  Nevertheless  still  another  proof  was  carried  out. 

It  was  based  on  the  fact  that  the  scattering  ratios  considered  as  a  function 
of  q  should  extrapolate  convincingly,  at  q  *  O.to  the  data  reported  previously 
for  monodisperse  systems.  That  this  applies  in  fact,  is  illustrated  in 
Figures  3-5.  In  order  to  make  the  extrapolations  to  q  *  0  as  secure  as  possible 
the  interval  of  integration  had  to  be  reduced  from  to.  -  0.1  to  0.01. 

The  computation  of  the  £  values  Is  the  secct  -  step  in  the  overall 

computation  procedure,  the  first  step  being  the  evaluation  of  the  quantities 

Ifj  and  Ij__  given  In  Equations  (3a)  and  (3b).  Essentially,  all  statements  made 

with  regard  to  the  computation  of  (T  apply  also  to  the  computation  of  these 

quantities.  Figure  C,  .hows  on  the  example  of  that  it  also  extrapolates  very 

well  to  q  =0.  The  values  pertinent  to  q  *  Ov  taken  from  the  hook  by 

4 

Denman,  .feller,  Pangonis  . 

17.  PRELIMINARY  DISCUSSION  OF  THE  RESULTS  OBTAINED. 

It  is  not  intended  here  to  discuss  the  results  proper, their 
characteristic  features  and  practical  application.  Ibis  will  be  done  In  a 
separate  article  to  be  published  elsewhere.  However,  a  few  salient  features 
may  he  referred  to  here.  First  of  all,  it  is  of  interest  to  consider  briefly 
the  character  of  the  6^ vs .  q  curves  given  in  Figures  3-5.  The  curves  jo J.ricri  w 
to  large  y  -values  (  y  -  90°  and  150°)  in  Figure  3  are  quite  simple;  this  is 
easy  to  understand  since  minima  and  maxima  of  S'  enter  the  forward  direction 
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/ 


O  1  p 

(7>90°)  only  at  rather  high  Ct  values.  The  re  lore,  an  increase  in 
beterodispersity  simply  increases  the  S -values  monotonically  at  these 

large  angles  of  observation.  Similarly,  the  curves  are  bound  to  be 
oonotonic  at  low  p  values  if  y <90°  as  long  as  q  is  reasonably  small,  the 
reason  being  again  the  absence  of  minima  and  maxima.  This  is  apparent  from 
Figure  k  v.  \h  applies  to  /  =  10° .  However,  at  larger  p-values  for  the 
came  y,  6  passes  through  a  minimum  or  a  maximum  within  the  same  range  of 
relatively  small  q -values .  This  follows  from  Figure  5.  It  is  finally,  of 
interest  to  point  out  the  most  characteristic  feature  of  the  angular 
scattering  functions  derived  for  heierodisperse  systems:  Figure  7  shows 
for  p  *  1.0  that  even  a  very  modest  degree  of  beteix--  l spersion  characterized 
hy  q  *  0.1  increases  the  numerical  value  of  6"  at  90°  by  more  than  The 

increase  is  about  100-fold  at  q  -  1.0.  In  addition,  there  is  an  8°  shift  In  the 
location  of  the  minimum  value  of  6  in  the  latter  instance,  It  follows  from  this 
comparison  that  the  analysis  of  the  size  distribution  curves  will  be  particularly 
attractive  at  those  angles  where  6" -minima  or  maxima  occur,  because  even  slight 
degrees  of  heterodlspersion  manifest  themselves  very  clearly. 

SUMMARY 

Fortran  programs  were  <1..  ,<-i  to  compute  light  scattering  functions  for 

beterodisperse  systems  which  can  be  described  by  a  bi-pa rametrlc,  unlmodal, 
exponential  distribution  function.  The  programs  had  to  include  the  evaluation 
of  the  Mie  functions  for  angular  scattering  for  any  arbitrary  values  of  the  size 
parameter  Cf  v.ithin  the  range  0  <  <2^25.  This  involved  the  generation  of  the 
Riccatl-Bessel  functions  and  Legendre  polynomials.  Two  methods  of  numerical 
integration  (viz.  Gaussian  Quadrature  and  Simpson's  rule)  of  the  expressions 
containing  the  products  of  the  Mie  functions  and  the  distribution  function  were 
investigated,  and  the  results  were  compared.  The  theoretical  results  were  aljo 
checked  by  extrapolation  to  the  condition  of  monodispersity  characterized  by 

q  — »0. 
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LEGENDS 


Figure  1 

SCHEMATIC  DIAGRAM  OF  THE  PROCEDURE  OF  COMPUTATIONS 

The  symbols  m,  y  and  a  of  the  text  are  replaced  here  by  M,  G  and 
A  due  to  technical  requirements  of  the  IBi  1403  printer  used.  Block  08, 
labelled  GROUP,  stands  for  a  group  of  subprograms  for  the  functions  Indicated. 

Figure  2 

VARIATION  OF  COMPUTED  ^-VALUES  WITH  THE  NUMBER  OF  NODES  ON 
USING  THE  GAUSSIAN  QUADRATURE  METHOD  ,  ... 

The  number  of  panels  Is  Indicated  as  the  parameter. 

m  »  1.20,  7  =■  90°,  p  =  q  »  4 .0 .  The  curve  was  traced  from  an  original  plotted 

by  the  Calcomp  digital  curve -plotter  with  a  separate  program  written  for  the 

purpose. 

3 

VARIATION  OF  6  WITH  q 
7  =  90°  and  150° 
m  =  1 . CO.  Parameter  p 

Figure  4 

VARIATION  OF  6  WITH  q.  I.  LOW  p -VALUES 
7  =  30° 

m  =  1.20. Parameter  p 


rcgure  l  m 

VARIATION  OF  E  WITH  q.  II.  INTERMEDIATE  p-VALUES 
Same  details  as  in  Figure  4. 


Figure  6 


VARIATION  OF  L,  WITH  q 
7  ~° 


m 


30^ 

1.20. Parameter  p 


Figure  7 

VARIATION  OF  (f  WITH  7  FOR  SMALL  DEGREES  OF  EETERODISPERSION 
{  ■  «  l.O)  'AND  p  =  1.0 


m  =  1.20 
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TABLE  1 

VALUES  COMPUTED  BY  (A)  GAUSSIAN  QUADRATURE  AND  (B)  SIMPSON'S  RULE 

m  -  1.20,  7  »  90 


q 

p 

*1 

*2 

B1 

B2 

h 

0.4 

i 

0 

.0000 

0.4 

4 

.31614 

.31613 

.3028 

2.0 

0 

.10107 

2.0 

4 

.41845 

.41847 

.41843 

.41846 

.4182 

4.0 

0 

.30938 

.31223 

.30720 

.30720 

.3049 

4.0 

4 

.58547 

.58547 

.57309 

.57309 

.5710 

*1  ! 
*2  S 
h  s 
B2  5 
V 


9  nodes 
5  nodes 

&x  ~  0.2 

££t  -  0.1 

Data  from  Ref.  8 


(computed  with 
Simpson's  rule, 
££t  =  0.2) 
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